Contact dynamics in a gently vibrated granular pile 
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We use multi-speckle diffusive wave spectroscopy (MSDWS) to probe the micron-scale dynamics 
of a granular pile submitted to discrete gentle taps. The typical time-scale between plastic events 
is found to increase dramatically with the number of applied taps. Furthermore, this microscopic 
dynamics weakly depends on the solid fraction of the sample. This process is strongly analogous 
to the aging phenomenon observed in thermal glassy systems. We propose a heuristic model where 
this slowing down mechanism is associated with a slow evolution of the distribution of the contact 
forces between particles. This model accounts for the main features of the observed dynamics. 
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Internal contact forces in dense granular systems are 
very inhomogeneous 0,0], even for crystalline assemblies 
3] . The stress at a given contact can change macroscop- 
ically following a relative displacement of the two parti- 
cles of order of microns. Hence, large modifications of 
the contact forces field can result from minute deforma- 
tions of the pile. This phenomenon is crucial in under- 
standing the catastrophic yielding occurring in granular 
systems submitted to a slowly varying stress (avalanches 
0, shear-bands in triaxial tests It also explains why 
the sound transmission through a granular sample can 
be strongly affected by very small deformations [6j. To 
probe the evolution of the internal stress field, one needs 
to measure forces directly 0] which is difficult in 3D 
. In this letter, we propose a different approach: we use 
MSDWS to measure particle displacements on micron- 
scales in a pile submitted to gentle discrete taps. These 
vibrations are too weak to induce large-scale rearrange- 
ments which would eventually lead to a compaction of 
the granular system [S3, UJJ, UT| • We can therefore evalu- 
ate the micro-dynamics of the contacts without significa- 
tively perturbing the packing structure. 

We use glass beads of diameter 45 ± 2 /im, contained 
in a glass cell (30 mm x 10 mm x 2 mm). To reduce 
electrostatic forces and the effects of moisture, the gran- 
ular system is saturated with pure water. During the 
experiment, the mean packing fraction <f> of is obtained 
by measuring the position of the upper surface of the 
pile with a CCD camera. Although a systematic error 
of 2% can not be avoided, this allows us to detect rela- 
tive changes in <j> as small as 0.01%. To produce motion 
in the pile, we use a piezoelectric actuator on which the 
cell is rigidly mounted. Vertical vibrations of precisely 
controlled amplitude, shape and durations can thus be 
applied to the granular column. In this experiment, we 
focus on a single type of mechanical excitation, later re- 
ferred to as a " tap" , which consists in a train of square 
wave vibrations of frequency 1 kHz and duration 100 ms. 
Different applied voltage are used, yielding various verti- 
cal amplitudes ranging from 50 to 300 run. 



In a standard experimental run, the pile is prepared by 
turning the cell upside down then allowing the particles 
to sediment for half an hour. This procedure yields re- 
producible structures of low volume fraction. The pile is 
then submitted to high amplitude taps (of vertical ampli- 
tude 300 nm) until it reaches a prescribed packing frac- 
tion cj) s . During this compaction stage, the evolution of 
the packing fraction 4> with the number of taps (Fig. 1) 
is highly reproducible, and consistent with previous ex- 
perimental results on dry granular systems. We however 
note that the packing fraction we reach is well above the 
close packing limit expected for a fully disordered pile 
(~ 0.64). This indicates that cristallisation does occur 
in our system under vibration. This first step is for us a 
mean to reproducibly prepare a granular sample of given 
packing fraction with essentially the same preparation 
history. We then start probing the dynamics of contacts 
by submitting the cell to very gentle taps (of amplitude 
50 nm). As shown in Fig. 1, this second step does not in- 
duce significative evolution of the packing fraction except 
for initially very loose packs. 




FIG. 1: Four compaction curves : (+) is only excited by gentle 
vibrations; The other systems are prepared with (x) 50, (□) 
150, (A) 8000 high impulsions pulses, before being excited by 
gentle vibrations. There is a systematic error of 2% on the 
measurements of the packing fraction. 
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To probe the microscopic dynamics induced by these 
gentle taps, we use MSDWS [H Q. This technique, 
which allows one to resolve sub-micron displacements, 
has been successfully applied to granular dynamics by 
several groups [ljj . The sample is illuminated with 
a He-Ne laser beam at a depth of 2 cm below the pack 
upper surface (1 cm over the bottom). Photons are mul- 
tiply scattered by the particles 0, and form a speckle 
pattern on the opposite cell wall that we record with a 
CCD camera. In the absence of vibrations, the speckle 
image does not change in time as temperature is insignif- 
icant for such large objects. By contrast, the taps induce 
some irreversible particles displacements which modify 
the speckle image. To quantify this internal dynamics, 
we measure the intensity correlation of speckle images, 
taken between taps, as a function of the number of taps t 
that separate them. This function generally depends on 
the total number of small amplitude taps t w that have 
been performed. We therefore calculate the two-times 
correlation fonction g(t w ,t): 



(I(t w + t) ■ I(t w )) sp kl 



spkl 



{I(tw) 2 ) spk i (I(tw))spkl 
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In this expression, ( ) sp ki corresponds to averaging over 
different speckles. MSDWS thus allows one to rapidly 
access relaxation time by substituting time- with space- 
averaging and is therefore well suited to the study of 
non-stationary dynamical systems. 

Figure [21 shows three correlation fonctions obtained 
with the same sandpile at different values t w . These fonc- 
tions are well fitted by stretched exponentials: 

g{t w ,t) = exp (- (t/r^)) ^ 




FIG. 2: Three different correlation fonctions obtained after 
three different times t w . Solid lines correspond to experimen- 
tal values and dashed lines to the stretched exponential fit. 

For different packing fractions, we follow the evolution 
of the dynamics by monitoring the two parameters r(t w ) 
and a(t w ) as a function of the total number of gentle 
pulses t w . We find that the exponent a is roughly 



constant (~ 0.8 ± 0.2) and independent of the packing 
fraction <p s . By contrast, the time r(t w ) increases by 
five decades over the range of t w explored, as shown in 
figure |3J It should be noted that this dynamics can be 
immediately reset by submitting the system to a few taps 
of larger intensity (such as those used for compacting 
the sample). A careful examination of the r(t w ) curve 
also reveals large fluctuations in the internal dynamics, 
especially in looser packs for which the packing fraction 
slowly evolves. During certain periods of time, the 
dynamics is restarted as shown by a sudden decrease 
of r(t w ). This may correspond to catastrophic failures 
of the pack structure, which compete against a global 
reinforcement of the granular contacts. 
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FIG. 3: Evolution of the dynamical time r with the number 
of low magnitude taps t w for different packing fractions <j> 3 



This result demonstrates that the response of a gran- 
ular system to small perturbations is strongly depen- 
dent on the history of its preparation (the number of 
applied taps t w ), and rather insensitive to the pack- 
ing fraction. More quantitatively, the typical relax- 
ation time associated with a given mechanical excitation 
evolves as a power law of the total number of taps t w 
'■ T (tw) ~ tw 2±0 ' 2 - This observation is strongly remi- 
niscent of the aging behavior recently observed in many 
glassy colloidal systems In these materials, 

the longest (a— Relaxation time is found to grow as a 
power law of the time since the system was left to rest 
after a rapid shearing. In spite of this strongly analo- 
gous behavior, the microscopic processes leading to this 
dynamical arrest are qualitatively different. In colloidal 
systems, stress relaxation occurs by thermally activated 
rearrangements of the structure. In granular materials, 
temperature is effectively zero and relaxation only results 
from externally applied vibrations which induce the slip- 
page of some contacts (the most fragile ones). 

The frequency of these slipping events is directly 
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probed by the MSDWS technique. Slipping events 
indeed result in some grain displacements (either trans- 
lations or rotations) around the broken contact, until a 
new equilibrated configuration is found. In this gentle 
vibrations regime, these events are too small to con- 
tribute in a significant way to the compaction process. 
This means that the associated grain displacements are 
of order 6 « D, where D = 45 /im is the particle di- 
ameter. Although tiny, these irreversible displacements 
control the speckle decorrelation. However, to get a 
quantitative estimate of their frequency, one needs to 
assume that they are uniformly distributed in space and 
have a unique characteristic amplitude 8. Within this 
hypothesis, r(t w ) is directly proportional to the inverse 
of the yielding frequency [12, bjj • 

We now turn to a tentative microscopic model to cap- 
ture this slowing down process. As underlined before, 
the observed dynamical arrest is rather insensitive to the 
packing fraction of the sample. We thus need to intro- 
duce another internal variable that would control the in- 
stantaneous response of the pack to gentle vibrations. 
Here we propose to focus on contact stress distribution. 
It has been observed that the form of this distribution 
is almost independent of the volume fraction of the pack 
and the preparation history. However, standard measure- 
ments are not sensitive enough to detect small variations 
in these distributions especially in the low force limit, 
that may follow from very gentle mechanical vibrations. 
We will argue here that the observed evolution of the 
dynamics results precisely from a slow modification of 
the stress distribution inside the pack which effectively 
rigidify the granular pile. 

To modelize such a dynamics, we picture the gran- 
ular assembly as a set of independent contacts (which 
total number is supposed to be a constant.) Each con- 
tact is characterized at time t by the normal and tan- 
gential component of the contact force which we denote 
a n and at respectively. Mechanical equilibrium imposes 
that a t < ^a n . For simplicity, we will make the friction 
coefficient /i equal to 1 in the rest of the letter. For a 
given packing structure, the state of the internal stress 
field is characterized by the two variable stress density 
distribution P <T (a n , a t ). As the cell is vibrated, mechan- 
ical waves travelling through the sample induce random 
stress fluctuations on each contact. Such perturbations 
can locally trigger the rupture of a contact, whenever the 
shear force a t overcomes the normal force a n . We assume 
an exponential distribution %(<5ct) of the maximum force 
fluctuation induced by the tap on each contact: 



x(Sa) 



1 



exp 



5a 



(2) 



In this expression, the mean force fluctuation Sa is an 
increasing function of the applied vibration amplitude. 



Thus the probability for a given contact to yield following 
a single tap writes: 



uj y (a n ,a t ) = exp == 

oa 



(3) 



This expression is a consequence of the peculiar form 
(Eq. [2J taken for the tap induced force fluctuations 
x(5a). However, the main results of the present model re- 
main valid for any fast decaying distribution (faster than 
a power law). 

After a yielding event, the force at the renewed con- 
tact is chosen from a given "rejuvenated" distribution 
which we consider intrinsic to the system. The distribu- 
tion of normal forces in a granular pile under moderate 
load is known to exhibit an exponential tail at high forces 
and a plateau below the mean force 0, 0, Q • Numerical 
measurements have also shown that, for a given value of 
the normal force, the tangential forces are uniformly dis- 
tributed between and the sliding limit a t = [ia n (in a 
2D case) [i^ • We use these different observations to infer 
the form of the "rejuvenated" distribution P re j(a ni a t ), 
which is thus written: 



P re j(a n , a t ) 



ao ■ a 7 , 



■ exp 



(4) 



where ao is the mean stress inside the pile. For simplicity, 
we have omitted the plateau saturation of the distribu- 
tion at low forces. We can now derive the dynamical 
equation of evolution of the stress distribution P a : 



dPq(a n ,at) 

at 



= -Pa{a n ,a t ) ■ Lo v (a n ,a t ) (5) 

L+ Prei{On,Ot)-F{P a ) 



where F(P a ) is the total frequency of sliding events which 
is self-consistcntly defined as: 



F{P°) = 



P a (a' n ,a' t ) ■ uj y (a' n ,a' t ) ■ da' n da' t (6) 



The present description exhibits many common fea- 
tures with Bouchaud's trap model [23| of glass transition. 
In the latter, the internal dynamics of a glassy liquid is 
pictured as a succession of thermal escapes from energy 
wells of various depths. In an analogous way, each con- 
tact here can be considered as frozen in a mechanical trap 
(the local solid friction cone), the depth of which depends 
on the relative amplitude of the normal and shear com- 
ponents of the contact force. Moreover, in the absence of 
temperature, mechanical vibrations plays the role of the 
energy source by allowing individual contacts to hop out 
of their trap. 

As in the trap model, we thus observe two limiting 
regimes depending on the relative values of the intensity 
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of the applied stress 5a and the width ag of the rejuve- 
nated distribution. For large vibrations, i.e. 5a > ao, 
the rejuvenated stress distribution is a stationary solu- 
tion of Eq. ©, and the yielding frequency is constant 
with time. By contrast, for 5a < ao, the stress distribu- 
tion keeps evolving endlessly. Figure 0] shows the time 
evolution of P a obtained by numerically solving Eq. (JBJ 
for da = co/20, starting with P a = P re j- It shows that 
fragile contacts - contacts of low normal force or close to 
the sliding limit (inset) - are slowly depleted. As a result, 
the number of sliding events per time unit decays. More 
quantitatively, we find that the characteristic time be- 
tween events grows linearly with the elapsed time. This 
is to be compared with the r ~ t^ 2 scaling behavior ob- 
served in the experiment (Figure [3]) 
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FIG. 4: Results of the numerical model for uq = 1 and 
Scr — cto/20: (a) Distributions of normal forces for t w ~ 
1,10,100,1000 and 10000 (from left to right); (b) distribu- 
tions of tangential forces for a n = ao at the same times t w 
(from bottom to top). 

We have evidenced, through MSDWS measurements, 
the existence of a slowing down of the micro-scale dy- 
namics over more than five decades in gently vibrated 
granular piles. This behavior is reminiscent of the aging 
process observed in glassy systems. This dynamics ap- 
pears to be weakly connected to the overall grain-scale 
structure, which suggests a two-level description of gran- 
ular systems. At high enough vibration, a granular pile 
evolves through the restructuration of the piling geom- 
etry, leading to a slow irreversible compaction. In this 
regime, forces networks are rapidly renewed and show 
no history-dependent behavior. At very low vibrations 
however, the geometry of the pile is essentially frozen, 
but the forces network can still evolve by slowly deplet- 
ing the most fragile contacts. This leads to an effective 
reinforcement of the pack structure as is evidenced in the 
present study by the decrease of vibration induced plastic 
events. 

The precise nature of the yielding events remains how- 
ever unclear. In particular, one might expect large spatial 
and temporal correlations between them, which we can- 
not probe with DWS. Another important question con- 



cerns the relevance of such modifications to the onset of 
macroscopic flow. For instance, does the observed rein- 
forcement of the force networks play a role in changing 
the threshold of avalanche triggering or shear-banding 
appearance ? 

We wish to thank Luca Cipelletti and Jean-Marc di 
Meglio for fruitful discussions and Deniz Gunes for help- 
ing us with the estimation of I*. 



[1 
[2. 
[3; 

[4; 

[5. 

[7; 

[8] 
[9] 

[10; 

[11 

[12. 

[13 

tw; 

[15 
[16 

[17; 

[is; 

[19 
[20 

[21 
[22 
[23 



Electronic address: georges.debregeas@college-de-france.fr 
D. M. Mueth, H. M. Jaeger and S. R. Nagel, Phys. Rev. 
E 57, 3164 (1998). 

G. L0voll, K. J. Mal0y and E. G. Flekk0y, Phys. Rev. E 
60, 5872 (1999). 

D. L. Blair N. W. Mueggenburg, A. H. Marshall, H. M. 
Jaeger and S. R. Nagel, Phys. Rev. E 63, 041304 (2001). 
L. Staron, J.-P. Vilotte and F. Radjai, Phys. Rev. Lett. 
89, 204302 (2002). 

J. Desrues, R. Chambon, M. Mokni and F. Mazerolle, 
Geotechnique 46, 529 (1996). 

C. -h. Liu, Phys. Rev. B 50, 782 (1994). 

D. Howell, R. P. Behringer and C. Veje, Phys. Rev. Lett. 
82, 5241 (1999). 

M. Da Silva and J. Rajchenbach, Nature 406, 708 (1999). 

J. B. Knight C. G. Fandrich, C. N. Lau, H. M. Jaeger 

and S. R. Nagel, Phys. Rev. E 51, 3957 (1995). 

P. Philippe and D. Bideau, Europhys. Lett. 60, 677 

(2002). 

C. Lesaffre, V. Mineau, D. Picart and H. Van Damme, 

C. R. Acad. Sci. Paris, t.l, Serie IV, 647 (2000). 
G. Maret and P. E. Wolf, Z. Phys. B. 65, 409 (1987). 

D. J. Pine, D. A. Weitz, P. M. Chaikin and E. Her- 
bolzheimer, Phys.Rev.Lett. 60, 1134 (1988). 

L. Cipelletti and D. A. Weitz, Rev. Scien. Inst. 70, 3214 
(1999). 

V. Viasnoff, F. Lequeux and D. Pine, Rev. Scien. Inst. 
73, 2336 (2002). 

N. Menon and D. J. Durian, Science 275, 1920 (1997). 
K. Kim, J. J. Park, J. K. Moon, H. K. Kim and H. K. Pak, 
Journal of the Korean Physical Society 40, 983 (2002). 
I* was found equal to 220 ± 20/^m, which allows in our 
geometry the use of a diffusive model for the light in the 
sample [l2||. 

A. Knaebel, M. Bellour, J.-P. Munch, V. Viasnoff, F. 
Lequeux and J. L. Harden, Europhys. Lett. 52, 73 (2000). 
V. Viasnoff and F. Lequeux, Phys. Rev. Lett. 89, 065701 
(2002). 

L. Cipelletti, L. Ramos, S. Manley, E. Pitard, D. A. 
Weitz, E. E. Pashkovski and M. Johansson, Faraday Dis- 
cuss. 123, 237 (2003). 

D. J. Durian, D. A. Weitz and D. Pine, Science 252, 686 
(1991). 

F. Radjai, M. Jean, J.-J. Moreau and S. Roux, Phys. 

Rev. Lett. 77, 274 (1996). 

J.-P. Bouchaud, J. Phys. I 2, 1705 (1992). 



